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Abstract. We establish a connection between the stabihty of an eigenvalue under a magnetic per- 
turbation and the number of zeros of the corresponding eigenfunction. Namely, we consider an eigen- 
function of discrete Laplacian on a graph and count the number of edges where the eigenfunction 
changes sign (has a "zero"). It is known that the n-th eigenfunction has n — 1 + s such zeros, where 
' the "nodal surplus" s is an integer between and the number of cycles on the graph. 

. We then perturb the Laplacian by a weak magnetic field and view the n-th eigenvalue as a function 

of the perturbation. It is shown that this function has a critical point at the zero field and that 
O , the Morse index of the critical point is equal to the nodal surplus s of the n-th eigenfunction of the 

' unperturbed graph. 

^ ■ 1. Introduction 

r-| , Studying zeros of eigenfunctions is a question with rich history. While experimental observations 
^ \ have been mentioned by Leonardo da Vinci [35], Galileo [23] and Hooke [9], and greatly systematized 
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by Chladni [M], the first mathematical result is probably due to Sturm [H]. The Oscillation Theorem 
of Sturm states that the number of internal zeros of the n-th eigenfunction of a Sturm-Liouville 
operator on an interval is equal to n — 1. Equivalently, the zero of the n-th eigenfunction divide 
I> ■ the interval into n parts. In higher dimensions, the latter equality becomes a one-sided inequality: 
Courant [TTf 118] proved that the zero curves (surfaces) of the n-th eigenfunction of the Laplacian 
divide the domain into at most n parts (called the "nodal domains"). 
^ ] Recently, there has been a resurgence of interest in counting the nodal domains of eigenfunctions, 
O [ with many exciting conjectures and rigorous results. The nodal count seems to have universal features 
[121 [131 EZ], is conjectured to resolve isospectrality [21], and has connections to minimal partitions 
of the domain [23 [7], to name but a few. For a selection of research articles and historical reviews, 
see the dedicated volume 

^ ; On graphs, the question can be formulated regarding the signs of the eigenfunctions of the operator 

c5;(l) if ->MI^I, H = Q-C, 

where V is the set of the vertices of the graph, Q is an arbitrary real diagonal matrix and C is 
the adjacency matrix of the graph. On a graph, by a "zero" we understand an edge on which the 
eigenfunction changes sign, and not the exceptional (with respect to perturbation of Q) situation of 
an eigenfunction having a zero entry. 

The subject of sign changes and nodal domains (connected components of the graph left after 
cutting the above edges) was addressed in, among other sources, Fiedler ^2] who showed the analogue 
of Sturm equality for tree graphs (see also [lOj), Davies, Gladwell, Leydold and Stadler [19j, who 
proved an analogue of Courant (upper) bound for the number of nodal domains, Berkolaiko [6], who 
proved a lower bound for graphs with cycles and Oren [38] , who found a bound for the nodal domains 
in terms of the chromatic number of the graph. A number of predictions regarding the nodal count 
in regular graphs (assuming an adaptation of the random wave model) is put forward in [20]. For 
more information, the interested reader is referred to the book [11] and the review [3]. 
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The study of magnetic Schrodinger operator on graphs has a similarly rich history. To give a sam- 
ple, Harper [25] used the tight-binding model (discrete Laplacian) to describe the effect of the mag- 
netic field on conduction (see also [30]). In mathematical literature, discrete magnetic Schrodinger 
operator was introduced by Lieb and Loss [M] and Sunada [13,116], and studied, among other sources, 
in [m US] (see also [17] for a review) . 

In this paper we present a surprising connection between the two topics, namely the number of sign 
changes of n-th eigenfunction and the behavior of the eigenvalue A„ under the perturbation of the 
operator if by a magnetic field. To make a precise statement, we need to introduce some notation. 

Consider a generic eigenfunction on the graph, that is an eigenfunction that corresponds to a 
simple eigenvalue and is nonzero at the verticeso We denote by the number of sign changes 
(also called sign flips, hence the notation 0) which are defined as the edges of the graph at whose 
endpoints the eigenfunction has different signs. Here n is the number of the eigenfunction in the 
sequence ordered according to increasing eigenvalue. The combined results of [221 M E] bound the 
number by 



where /3:=|V^| — |i?| + lis the first Betti number (number of independent cycles) of the graph. Here 
and throughout the manuscript we assume that the graph is connected. We will call the quantity 



the nodal surplus. This is the extra number of sign changes that an eigenfunction has due to the 
graph's non-trivial topology. 

Magnetic field on discrete graphs has been introduced in, among other sources, [311 113 US]. Up 
to unitary equivalence, it can be specified using /3 phases a = (aj)^^-^ G (— vr, vr]^. We consider the 
eigenvalues of the graph as functions of the parameters a. The zero phases, a = 0, correspond to 
the graph F without the magnetic field. We are now ready to formulate our main result. 

Theorem 1.1. The point a = is the critical point of the function A„(a). // A„(0), the n-th 
eigenvalue of the non-magnetic operator on T, is simple and the corresponding eigenfunction has no 
zero entries, its nodal surplus an is equal to the Morse index — the number of negative eigenvalues 
of the Hessian — of A„(a) at the critical point a = 0. 

An immediate consequence of this theorem is the following. 

Corollary 1.2. The non- degenerate n-the eigenvalue of the discrete Schrodinger operator is stable 
with respect to magnetic perturbation of the operator if and only if the corresponding eigenfunction 
has exactly n — 1 sign changes. 

By "non- degenerate" we understand a simple eigenvalue whose eigenfunction does not vanish on 
vertices and by "stability" we mean that the eigenvalue has a local minimum at zero magnetic field. 

Other possible consequences of our result and links to several other questions are discussed in 
Section The rest of the paper is structured as follows. In Section [2] we provide detailed definitions. 
Section [3] is devoted to a duality between the magnetic perturbation and a certain perturbation to 
the potential, coupled with removal of edges. This leads to an alternative proof of the result in the 
case /3 = 1 (subsection 13.31) which, although unnecessary for the general proof, provides us with 
some important insights. Section HI collects the tools necessary for the proof of Theorem 11.11 while 
Section [5] contains the proof itself. 



This is the generic situation with respect to the perturbation of the potential Q. 



(2) 



n — l<0n<'ri — 1 + /3, 



(3) 
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2. Magnetic Hamiltonian on discrete graphs 

Let r = (V, E) be a simple finite graph with the vertex set V and the edge set E. We define the 
Schrodinger operator with the potential q : V M.hj 

(4) :MI^I ^RI^I, (H^)^ = -J2^v + qu^u, 
that is the matrix H is 

(5) H = Q-C, 

where Q is the diagonal matrix of site potentials qu and C is the adjacency matrix of the graph. It is 
perhaps more usual (and physically motivated) to represent the Hamiltonian as H = Q + L, where 
the Laplacian L is given hj L = D — C with D being the diagonal matrix of vertex degrees. But 
since we will not be imposing any restrictions on the potential Q, we absorb the matrix D into Q. 
The operator H has \V\ eigenvalues, which we number in increasing order, 

Ai < A2 < ... < A|y|. 

We define the magnetic Hamiltonian (magnetic Schrodinger operator) on discrete graphs as 

(6) (i/V')« = -5Ze*^-V^„ + g„7/>„, 

with the convention that A^^^ = ~A.u,v, which makes H self-adjoint. For further details, the reader 
should consult [Ml SHI [151 • 

A sequence of directed edges C = [ui,U2, ■ ■ ■ , Un] is called a cycle if the terminus of edge Uj coincides 
with the origin of the edge Wj+i for all j (u„+i is understood as mi). The fiux through the cycle C by 

(7) $c = ^ui,u2 + ■■■ + ^«„-i,«„ + Au^,ui mod 271 

Two operators which have the same fiux through every cycle C are unitarily equivalent (by a gauge 
transformation). Therefore, the effect of the magnetic field on the spectrum is fully determined by 
/3 fluxes through a chosen set of basis cycles of the cycle space. We denote them by ai, . . . , a/3 and 
consider the n-th eigenvalue of the graph as a function of a. 

More precisely, flx an arbitrary spanning tree of the graph and let 5* be the set of edges that do 
not belong to the chosen tree. Obviously, S contains exactly (3 edges. 

Lemma 2.1. Any magnetic Schrodinger operator on the graph T is unitarily equivalent to one of 
the operators of the type 

u = v, 

(8) H^,,= {-1, {u,v)eS\S, 

-e^*°% {u,v) = seS, 

where the sign of the phase is plus if u < v and minus if u> v. 

Example 2.2. Consider the triangle graph — a graph with three vertices and tree edges connecting 
them. One of the equivalent forms of the magnetic Hamiltonian for this graph is 

j qi -e^" -1' 
^(r^g) = -e"^" q2 -1 
V -1 -1 93 

The spectrum of HiT'^ ) as a function of a E (— 7r,7r] is shown in Fig. 
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Figure 1. The eigenvalues of the triangle graph as functions of a magnetic phase a 
(bold lines) and the eigenvalues of the unperturbed graph (horizontal lines). 

3. A DUALITY BETWEEN A MAGNETIC PHASE AND A CUT 

In this section we explore a simple result which shows a connection between two types of pertur- 
bations of the operator H that will be used to prove the main theorem. It illustrates the duality 
between the perturbation of a discrete Schrodinger operator by a magnetic phase on a cycle and the 
operation of removing ("cutting") an edge that lies on the cycle. The latter operation was used to 
prove the lower bound on the number of nodal domains in [6j and to study partitions on discrete 
graphs in [8]. 

3.1. Tools used. The result of this section (Theorem 13.31 below) is based on the following version 
of the Weyl's inequality of linear algebra that can be obtained using the variational characterization 
of the eigenvalues (see [211 Chap. 4] for similar results). 

Theorem 3.1. Let A he a self-adjoint matrix and B be a rank-one positive semidefinite self-adjoint 
matrix. Then 

(9) \n{A -B)< K{A) < Xn+i{A - B), 

where A„ is the n-th eigenvalue, numbered in increasing order, of the corresponding matrix. Moreover, 
the inequalities are strict if and only if Xn{A) is simple and its eigenvector is not in the null-space of 
B. 
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Similarly, when B is negative-definite, we have 

(10) \n-i{A-B)<\n{A)<\niA-B), 

with an analogous condition for strict inequalities. 

Another useful result is the first term in the perturbation expansion of a parameter-dependent 
eigenvalue. Let A{x) be a Hermitian matrix-valued analytic function of x. Let A(a;) be an eigenvalue 
of the matrix A that is simple in a neighborhood of a point xq. We know from standard pertur- 
bation theory [33] that X{x) is an analytic function. Denote by u{x) the normalized eigenvector of 
corresponding to the eigenvalue A and by v all other normalized eigenvectors (in a slight abuse of 
notation). Then we have the following formula for the derivative of A evaluated at the point x = Xq. 

(11) ±x=/u—u 

dx \ ' dx 

3.2. Two operations on a graph. Let A„ be a simple eigenvalue and the corresponding eigenfunc- 
tion / be non-zero on vertices. Let (mi, U2) be an edge that belongs to one of cycles of the graph. We 
allow the graph to have magnetic phases on some edges, but assume that there is no phase on the 
edge {ui,U2). Then the operator H = Q — C has the following subblock corresponding to vertices Ui 
and U2, 



(12) i^(r)[„„ 



-1 

-1 



We consider two modifications of the original graph. The first modification of the graph is a cut: 
we remove the edge (ui,M2) and change the potential at sites ui and U2- Namely, we change the 
[ui, U2] subblock to 



(13) H{T 



cut\ 

7 )[ui,U2] 



9«i - 7 
- 1/7 



and leave the rest of the matrix H intact. We denote this modification by if(r™*). Note that this 
modification is a rank-one perturbation of the original operator H(r). Let B be the perturbation 
such that if(r^"*) = H{r) — B'^. Namely, the matrix B'^ has the [■Ui,'U2] subblock 



(14) b: 



Ul,U2\ 



7 -1 

-1 1/7 



and the rest of the elements are zero. Then B'^ is positive-definite if 7 > and negative-definite if 
7 < 0. Note that the cases 7 = 00 and 7 = can also be given meaning of removing (or imposing the 
Dirichlet condition at) the vertex ui or the vertex U2 correspondingly. However, we will not dwell on 
this issue and exclude these cases from our consideration. 

Notably, if / is an eigenfunction of H(T) and 7 = /uj/Zui ^ K, then / is also an eigenfunction of 
i/(r^^"*). Equivalently, / is in the null-space of the perturbation 5^. 

The second modification of the original graph is the introduction of a magnetic phase on the edge 
{ui,U2). The ['Ui,M2] subblock of the new operator H(T'^^^) is 

(15) -^(rmag)[Ml,M2] ~ ^ — e"**^ q 

while other entries coincide with those of H{r). Note that if(r°^g) is not a rank-one perturbation 
of H{r). However, it is a rank-one perturbation of the cut graph H{r™^) for any values of a and 7. 
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Namely, i/(r^"*) = //(F" ) - 5™^ where 



(16) B 



mag 

7 



mc 

ui,U2] ~ \ _p-iaj ]_^^ 



and all other entries of B'^^ are zero. Also, the spectrum of HiV^^^) and H(T) coincide when a = 
since the operators coincide. 

3.3. A duality between the two operations. We now want to apply Theorem 13. II to the spectra 
of r, r™* and r°j^g. However, we must take care to distinguish the two cases that correspond to 
equations ([9]) and (fTOj) (7 > and 7 < correspondingly). 

Definition 3.2. The eigenvalues of F, F™* and F^^^g will be numbered in increasing order starting 
from 1. When we happen to index "nonexistent" eigenvalues we use the following convention: 

A,(F) : 



-00, J < 1, 

00, j > n. 



Theorem 3.3. Let p{'~f) be 1 if < and otherwise. Then the following inequalities hold 

(17) A„_p(^)(F^^''*) < A„(F^3^g) < A„_p(^)+i(F!^"*), 
for all values of a and 7. Furthermore, for any fixed n 

(18) max A„_p(^) (F^*) = mm A„(F^,g) 

and 



cut\ 



(19) max A„(F°^g) = mm A„_p(^)+i(F: 

Finally, if there are no magnetic phases on the graph F (i.e. all entries of H{T) are real), then one of 
the extrema ( fTSl) or ( !T9l) is equal to A„(F) = A„(F^'^g), while the other is equal to A„(r) := A„(F"^'^g). 



Remark 3.4. Note that at this point we don't know which extremum, (JT8l) or ( JT9l) . is equal to A„(F). 
This information is related to the nodal surplus. We have also defined yet another modification of 
the graph F, the graph F whose adjacency matrix has —1 in place of 1 for the entries C^j^^j 

Cui ,U2 • 

Remark 3.5. Let M = MU{— 00, 00} be the extended real line and R = ]R/[— 00 = 00] be its projective 
("wrapped") version. The eigenvalue A„_p(^)(F™*) is then a continuous function of 7, considered as 

a function from M to M. See Figure |2] for an example. Note that according to our definitions, 
An-p(7)(FJ^"') = —00 for 72 = 1 and 7 < 0. 



Proof of Theorem \3.3[ The inequalities follow directly from Theorem 13. H since for any a the graph 
F°jjg is a rank-one perturbation of V^^. Whether it is positive- or negative-definite depends on the 
sign of 7, and results in the shift by p. 

The properties of the extrema we can get as follows. First of all, observe that if max An-_p(F™ ) — 
min A„_p+i(FJ^'^*), then XniXmag) is constant and equal to the common value of A„_p(F;^"*) and 

Let now max Xn—p{T^^) < min Xn-p-\-i{T'^^) . The eigenvalues of a one-parameter family can always 
be represented as a set of analytic functions (that can intersect). Let A'(F™*) be the analytic function 



NODAL COUNT OF GRAPH EIGENFUNCTIONS VIA MAGNETIC PERTURBATION 



7 




Figure 2. The duality between a magnetic field on one side and cut edge with added 
potential on the other. The graph is a triangle. The black curves correspond to the 
eigenvalues as functions of the magnetic phase. Colored symbols correspond to varying 
potential after cutting the edge. The x-axis ranges from — 7r/2 to 7r/2 with the magnetic 
phase taken as a = 2x and the potential parameter 7 = tan(x). The horizontal blue 
lines are eigenvalues of the true graph, while horizontal red lines are eigenvalues with 
magnetic phase vr. 

that achieves the maximum maxA„_p(r™*) and / be the corresponding eigenfunction. We will 
differentiate A'(rj^"*) using equation ( ITTl) . At the maximum point 7 = 7 we have (see equation ( fT3l) ). 



(20) 

From here it follows that 
(21) 7: 

Let 5 be the solution of 







dX' 
(^7 



07 



U2 I 



l/wi I 
iful I fu 



or, equivalently, 1 7/^, //^^ | = 1 . 

. Direct calculation shows that the eigenfunction / is in the 



null-space of the perturbation B"^'^ of f ll6l) with a = a and therefore / is both in the spectrum of F 



cut 



mag 5 



implying equation ( IT8|) follows. Proof of equation ( !T9|) is completely 

leading to the condition 



and in the spectrum of F 
analogous. 

Note that we could instead differentiate the eigenvalue of F^^^g, 



(22) 



f f e'°' G 
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instead of equation ( l20l) . One then sets 7 = e^'^fu^lfu^ G M to the same effect. 

Finally, when the matrix if(r) is real, the eigenfunctions of F™*, F^'^g and F^'^g are real-valued. 
When a = we can verify directly that the eigenfunction / of FJJ^^q is also an eigenfunction of F~"* 
by setting 7 = juil Jux- When a = vr, we also set 7 = fuilfux ^-iid do the same. □ 

Theorem 13.31 highlight a sort of duality between the two modifications of the graph F. The spectra 
of the graphs with a magnetic phase form bands (as the phase is varied) while the spectra of the 
graphs with the cut fill the gaps between this bands. Minimums of one correspond to maximums of 
the other and in half of the cases correspond to eigenvalues of the original graph. 

We now explain how the /3 = 1 case of Theorem 11.11 follows from Theorem 13.31 While for general 
/9 the proof is significantly different (it bypasses the interlacing inequalities and goes straight to the 
quadratic form), some key features are the same as in this simple case. 

Starting with the eigenvalue A„ of F and the corresponding eigenfunction /, we cut an edge on the 
only cycle of F to obtain a family of trees F™*. For 7 = 7 = /Ma/Ai, we have either 

maxA„_p(^)(F7*) = A„_p(^)(Ff ) = A„(F) = minA„(F^,g), 

or 

maxA„(F" ) = A„(F) = A„_p(;^)+i(Ff ) = min A„_p(^)+i(F7*). 

In the first case, according to Fiedler theorem (equation ([2]) with /3 = 0), the function / has ri— ^(7) — 1 
sign changes with respect to the tree F™*. Adding back the removed edge [ux.U'i) adds another sign 
change if 7 < and doesn't change the number of sign changes otherwise. In other words, it adds 
p(7) sign changes. Thus, with respect to F, the function / has n — \ sign change and (j„ = 0. In the 
second case, we similarly conclude that / has n — p{^) sign changes with respect to and n sign 
changes with respect to F. The nodal surplus is (T„ = 1. 

On the other hand, in the first case A„(F) is a minimum of A„(F|^^g) (Morse index 0), while in the 
second it is a maximum of A„(F^j^g) (Morse index 1), which shows that the Morse index coincides 
with an in the case (3 = 1. 

4. Tools of the main proof 
In this section we collect some basic facts that will be repeatedly used in the proof of Theorem ll.il 
4.1. Critical points of the quadratic form. 

Definition 4.1. Let F : M"' — )■ M be a twice differentiable function. If c is a critical point (i.e. 
V/(c) = 0), the inertia of c is the triple (n_,no,Ti+) that counts the number of negative, zero and 
positive eigenvalues correspondingly of the Hessian (the matrix of second derivatives) at the point c. 
The number n_ is called the Morse index (or simply index). 

The next theorem is a reminder that the eigenvectors of a self-adjoint matrix are critical points of 
the quadratic form on the unit sphere. 

Theorem 4.2. Let A be andxd real symmetric matrix and h{x) = {x,Ax), x G M*^, be the associated 
quadratic form. Then the (real) eigenvectors of the matrix A are critical points of the function h{x) 
on the unit sphere \\x\\ = 1. 

Let A„ be the n-th eigenvalue of A and Z^"-* be the corresponding normalized eigenfunction. Define 

(23) = ij-{Xm < A„}, = ij-{Xra = A„, n], n+ = ij-{Xm > A„}, 

with n^+no+n+ = d—l. Then the inertia of the critical point x = f^'^^ is (n_,no,n+). In particular, 
if A„ is a simple eigenvalue, the inertia is {n — l,0,d — n). 
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Remark 4.3. The value of the quadratic form h at the critical point /'•"^ is A„. 

Proof. The idea is intuitively clear: n__ — which is the Morse index — counts the number of directions 
in which the quadratic form decreases relative to the value at a; = These directions are the 

eigenvectors corresponding to the eigenvalues that are less than A„. Similar characterizations are 
valid for no and n_|_. 

We note that by Sylvester law of inertia, the inertia is invariant under the change of variables. 
Making the orthogonal change of coordinates to the eigen-basis of the matrix A, the quadratic form 
h{a) becomes 

h{a) = Xial + X2al + . . . + A^a^, 
while the sphere is given by the equations 

al + al + . . . + = 1. 

Thus, on the sphere, the quadratic form in terms of variables ai, . . . , a„_i, a„+i, . . . , is given by 

h = ^(Aj - A„)aJ, 

and the Hessian is a diagonal matrix with Xj — Xn,j = ^, ■ ■ ■ ,d, j n. The statement of the Theorem 
follows immediately. □ 

In the more general case when the matrix is Hermitian we should consider the quadratic form on 
the space C^. However, we can also consider it on the real space of double dimension. 

Theorem 4.4. Let A be an d x d Hermitian matrix and h{z) = {z,Az), z G C^, he the associated 
quadratic form. Consider h as a function of 2d real variables {x,y), where x G M*^, y E M.'^ and 
z = X + iy. Then the eigenvectors of the matrix A are critical points of the function h{x + iy) on the 
unit sphere \\x + iy\\ = 1. 

Let A„ be the n-th eigenvalue of A and f^"^ be a corresponding eigenf unction. Then the inertia of 
the critical point x+iy = /'■"^ with respect to the real space M^*^ is (2n_, 2no+l, 2n_|_), where n_, Uq and 
are defined by (!23|) . In particular, if Xn is a simple eigenvalue, the inertia is {2n — 2, l,2d — 2n). 

Proof. We adopt the convention that in the scalar product {z, w) the conjugation is applied to the 
first argument. Also, when mapping a 2; G to M^*^, we first list the real parts of the components 
of z and then the imaginary parts. With this conventions, it's an easy exercise to show that the 
quadratic form h{z) in variables (x, y) corresponds to the matrix 

(ReA -lmA\ 
^24) ^-[ir^A Re A J' 

Note that because A is Hermitian, the matrix Im A is scew-symmetric and, therefore, the matrix B is 
real symmetric. Every eigenvector z of the matrix A corresponds to 2 real eigenvectors of B with the 
same eigenvalue, namely (Rez,lmz) and {Iie{iz),lm{iz)) = (— Im^;, Re^;). These eigenvectors are 
orthogonal to each other, and to other similarly obtained eigenvectors. We therefore conclude that 
the spectrum of B is the same as the spectrum of A with all multiplicities doubled. The statement 
now follows from Theorem 14.21 □ 

4.2. Reduction to critical manifold. The tool introduced in this section is a simple idea already 
used in [21 [TJ E]. If we have a function /(xi, . . . , x„) with a critical point c then under some general 
conditions there is a (n — 1) -dimensional manifold around the point c on which the local minimum 
of / is achieved when we vary the variable Xi and keep the others fixed. Then the Morse index of / 
restricted to this manifold is the same as the Morse index of the unrestricted function. On the other 
hand, if the manifold is the locus of local maxima with respect to the variable xi, the Morse index 



10 



G. BERKOLAIKO 



on the manifold is one less than the unrestricted Morse index. The following theorem is a simple 
generalization of this idea. The proof is adapted from [7] . 

Theorem 4.5 (Reduction Theorem). Let X = Y (B Y-^. Let f : X ^ M. be a smooth functional such 
that (0,0) & X is its non- degenerate critical point with inertia Xx- Further, let, for every y &Y 
locally around 0, the functional f{y,y') as a function of y' has a critical point at y' = with inertia 
Xy±, that (locally) does not depend on y. Then the Hessian of f is reduced by the decomposition 
X = Y (B Y^ and the inertia of f with respect to the space Y is 

(25) Xy=Xx-Xyx. 

Proof. We calculate the mixed derivative of / with respect to one variable from Y and the other 
from Y^. In a slight abuse of notation we denote these variables simply by y and y'. We have 



^''^ dy 



df 



= 0, 

y=o 

since = is the critical point of f{y, y') as a function of y' for every y. Thus the Hessian of / has 
a block- diagonal form with two blocks that correspond to Y and Y-^. The spectrum of the Hessian 
is the union of the spectra of the blocks and the inertia is the sum of the inertias of the blocks, 

Ix = Xy + XyX . 

Equation (125|) follows immediately. □ 

The restriction of the function / to the subspace Y will be called the reduction of / with inertia 
Xy±. More often than not, we will be concerned with only the index of a critical point (the entry 
n_ of the inertia). In this case we will say "reduction with index m". The space Y is the locus of 
critical points of / with respect to the variable y' and we will refer to it as a critical manifold. 

Remark 4.6. Theorem 14.51 can be simply extended to the case when, for every fixed y, the critical 
point with respect to y' is located at y' = q{y) (rather than y' = 0). The function q{y) defines the 
critical manifold Q = {y, q{y)). If q{y) is a smooth function of y and g(0) = 0, the change of variables 

y ^y, y' h^y' - q{y) 

is non-degenerate (its Jacobian is a triangular matrix with Is on the diagonal) and makes / satisfy the 
assumptions of Theorem 14.51 By Sylvester law of inertia, the conclusion of the Theorem is invariant 
under the change of variables, that is inertia at point of / 1 ^ is 

= — Xy± . 
5. Proof of the main theorem 

We prove the main result in three steps. First we show by an explicit computation that the point 
is the critical point of the function A„(q;), where a = (ai, . . . ,0/3) G {—Tr,Tc]'^ are the magnetic 
phases. 

Then we fix an eigenpair A = A„(r) and /. We cut /3 edges of the graph turning it into a tree T, 
but modifying the potentials so that the eigenfunction / is also an eigenfunction of the tree T. It 
now corresponds to an eigenvalue number m, that is XmiT) = A. Considering the eigenvalue XmiT) 
as a function of the potentials, we find its inertia. This is done by considering the inertia of the 
corresponding quadratic form in the real space. The result of this step is related to the results on 
critical partitions, [8]. 

Finally, we relate the inertia of the function Xm{T) to the inertia of the function A„(a) at the 
corresponding critical points. This is done by relating the inertias of the quadratic forms, but now 
in the complex space represented as a real space of double dimension. 
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We recall that 5 is a set of f3 edges whose removal turns the graph T into a tree. By rmag'""'' 
we denote the graph obtained from T by introducing magnetic phases ai, . . . , on the edges from 
the set S. Similarly, by TJ^"* .^^ we denote the tree graph obtained by cutting every edge from S in 
a manner described in section 13.21 (see equation ( [T3l) and around). For future reference we list the 



7iv,7/3 

form that highlights the similarities between the three. 



quadratic forms of the original graph, the graph TJ^"* and the graph Fmag ' • We put them in the 



(27) 


h{x) 




^ ^ 2x^x^ 








u 


(n,i))e£\S 


{u,v)es 


(28) 














u 


{u,v)eE\S 


ej = {u,v)£S 


(29) 


h'''^^ (z) 




- E 2Re 


(^^^) - Y 2Re (z;:e'"^z„) 






u 


{u,v)£E\S 


ej=(?i,?))GS 



It is essential that the last quadratic form be considered on the complex vector space (that can 
be identified with a real space of double dimension). The first two forms can be considered on both 
real and complex spaces, with obvious modifications in the latter case. 

5.1. Critical points. Let / be an eigenfunction of the graph F. We have seen in Theorem 13.31 and 
its proof that the points a = and 7 = 7 (see equation fl2T|) ) are special: at these points / is an 
eigenfunction of the graphs T'^^^ and F™*. Moreover, they are critical points of the corresponding 
eigenvalues considered as function of the parameters a and 7 respectively. The result of this section 
generalizes this observation. 

Theorem 5.1. Let f be an eigenfunction of H{T) that corresponds to a simple eigenvalue A = A„(F). 
Assume f is non-zero on vertices of the graph F. For every edge {uj, Vj) ES,j = l,...,P, let 

(30) 7i = A//n,- 

Let p denote the number of negatives among the values 7^, 

P = Mlj<0, j = l,...,/3}. 

Then 

(31) A„(F) = A<^_,+i {^t.:'r,^ 

where (l)n is the number of sign changes of f with respect to the graph F. Moreover, the point 
(71, . . . ,7^) is a critical point of the function X^^-p+i (^r^^^* . 
Similarly for Fmag"'"'', 

(32) A„(F) = A. (F^--°) 
and (0, . . . , 0) is a critical point of the function A„(Fmag''°''')- 

Proof. It can be verified directly that / is an eigenfunction of the graph F-"* The nodal bound 
(j2]) with /3 = (proved by Fiedler in [22], see also [6]) shows that the eigenvalue corresponding to 
the function / has number /i' + 1 in the spectrum of the tree F-'^'* . Here /i' is the number of sign 
changes of / with respect to the tree. In general, this number is different from 0„ because we might 
have cut some of the edges on which / was changing sign. However, according to ( 130|) . these edges 
gave rise to negative values of 7^, therefore /x' = 0„ — p, proving equation (I3T]) . Equation (!32|) is 
trivial since T^--'^ = F. 

mag 
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To prove criticality of the points, we calculate the derivatives. Because the corresponding eigen- 
value is simple (for the tree, this follows from a theorem of Fiedler [22] and the fact that / does not 
vanish on vertices), the corresponding eigenvalues are analytic functions of the parameters and can 
be differentiated. 

Derivative of A(^^_p+i ^r™* j with respect to 7^ has been calculated in equation (120|) . resulting 



m 



(33) — A^„_p+i [T"^^ 

' (7iv,7/3) 

where we used the definition of 7, equation ( !30ll . 

The derivative of (rmag '"^) can be evaluated similarly using ( ITTj) . leading to 



7iv,7/3 



-i/.,r+i/.,r/7i = o. 



(34) Aa„ 



(o,...,o) 



since the eigenfunction / is real-valued. □ 

5.2. Index of the eigenvalue on the tree. In this section we elaborate on the first part of the 
result of Theorem 15. namely that (71, . . . ,7^) is a critical point of the function X^„^i (^r™* 

Theorem 5.2. Let f be an eigenfunction of H(T) that corresponds to a simple eigenvalue A = A„(r). 
Assume f is non-zero on vertices of the graph T and has 0„ sign changes. For every edge {uj, Vj) G S , 
J = 1, . . . , /3, let 

(35) lj = fvlfu- 

As before, p denotes the number of negatives among the values 7j. Then the point (71, . . . ,7^) as a 
critical point of the function A^^-p+i (^r^j,"* j has index n — 1 + (3 — (f)n- 

Proof. Denote by d the number of vertices of the graph F. Consider h'^^ (x), which is the quadratic 
form on the Hamiltonian of F™* as a function of d + P real variables {xi, . . . , Xd, 71, ... , 7/3) on 
the manifold x\ + . . . + x^ = 1. We note that the point (/i, . . . , Z^, 71, . . . ,7/3) is a critical point 
of h"^^ 7^3(^)5 can be easily shown by explicit computation. Indeed, the value of the Lagrange 
multiplier is the eigenvalue A„ and the gradient of 

F(a;i,...,Xrf,7i,...7^) = /i^"*..,^^(x) - A„(xi + . . . + x^) 

is zero: the first d equations become the eigenvalue condition Hf = Xnf and the last /3 are the same 

as das]). 

Denote the index of the point (/i, . . . , /rf, 71, . . . , 7^) by M. For every value of (71, ... , 7^) locally 
around the point (71, . . . ,7/3) the (0„ — p + l)-th eigenvector /^J^*.. ,^^ of F^j.'^'* is a critical point of 
^71* 7/3 ^ function of x. According to Theorem \4.2\ the critical point has index 0„ — p. 

According to standard perturbation theory (see, e.g., [33]) the eigenvector /^"* .^^ is a smooth 
(indeed, analytic) function of (71,..., 7^). This allows us to use Theorem 14.51 via Remark 14. 6 [ 
concluding that the critical point (71, . . . ,7^) of ^Jyi*..,^^(/^"*.._^^) has index M — (0„ — p). At this 
point we observe (see Remark 14. 3p that 



which is the function whose index we strive to evaluate. 



cut 

71. ■■■,7/3 
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eigenvector 
minimax, 

Figure 3. Schematic diagram of the proof of Theorem 15. 2[ The reductions are indi- 
cated by arrows, with the description of the parameters that are being reduced and the 
index of the reduction. Since we know the index of the critical point of h{x), we can 
follow the diagram, applying the Reduction Theorem, to calculate the index of 

\ / -pcut \ 

^<t>r,-p+l I J- 71, ...,7/3 )■ 



mm max 




On the other hand, consider x varying locally around the point /, so that the elements of x remain 
bounded away from zero. For each fixed x we look for a critical point with respect to the variables 
(71, ... , 7^). The terms of /i™* 7^(^) ^^^^ depend on a given 7 have the form 

(36) T(7) = -7x^ - x^/7. 

The critical point is (? = g{x) = x^/xu, which is a smooth function of x. The points (xi, . . . , Xd, gi, ■ ■ ■ 9/3) 
define another critical manifold to apply Theorem 14.51 to. Note that the critical manifold includes 
the point (/i, . . . , /^i, 71, . . . , 7^). Moreover, the critical point with respect to a given 7 is a maximum 
if (7 > and a minimum if g < 0. Each point is nondegenerate and, moreover, the sign of gj is locally 
the same as the sign of 7^ for all j. Different variables 7 are not coupled, thus the Hessian is diagonal. 
Therefore, the inertia of the points on the critical manifold is (/3 — p, 0,p) — it is a minimum with 
respect to p variables and maximum with respect to (3 — p. 

Consider now the function /i^y"* on the critical manifold. When 'j = g, the term (136|1 evaluates 

to 

T{g) = -2xuX,, 

and we find that, on the critical manifold, the function h™^ 7^(^) coincides with the quadratic form 
of the original graph, h{x). The point x = f, being the n-th eigenfunction of the graph is a critical 
point of h{x) and has index n — 1. Applying Theorem 14.51 we obtain 

n-l = M - {/3 -p). 

Coming back to (71, . . . ,7^3) as a critical point of A^^-p+i (^r™* we conclude that its index is 

M-{^n-p) = {n-l) + {f3-p)-{(Pn-p)=n-l + f3- 0,. 
The steps of the proof are summarized in Fig. 15.21 □ 

Remark 5.3. In [8] the eigenvalue of the tree graph F™* was interpreted as the energy of the 
"partition" with the given number of domains. Theorem 15.21 gives another route for the proof of the 
results of [8]. 

5.3. Index of the eigenvalue as a function of the magnetic field. Now we move from the 
critical point on the the tree to the critical point of the eigenvalue of the graph with magnetic 
phases. We apply the same method, practically retracing our steps, but now the quadratic form is a 
function of complex variables. 
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Figure 4. Schematic diagram of the proof of Theorem 15 .41 From Theorem 15.21 we 
know the index of the critical point of A<^„_p+i ^r™* j (bottom right corner). We 

then apply the Reduction Theorem four times to get the index of A„(rmag''"'') (bottom 
left corner). 

Theorem 5.4. Let f be an eigenfunction of H{T) that corresponds to a simple eigenvalue A = A„(r). 
Assume f is non-zero on vertices of the graph T and has (pn sign changes. Let Fmag"' ^ he the graphs 
with the magnetic phases ai,...,ap introduced on the edges from the set S. Then the index of 
(0, . . . , 0) as a critical point of the function A„(Fmag"' '') is the nodal surplus 0„ — (n — 1). 

Proof. Let z be the dimensional vector of complex numbers Zj = xj +iyj. We consider h'^^ ■y/si^)^ 
the quadratic form of the Hamiltonian of F™* , as a function of 2d + P real variables 

xi,...,Xd, yi,...,yd, 7i,...,7/3 

on the manifold + ... + l^;^^ = 1. We also consider hmag'°'^{z), the quadratic form of the 
Hamiltonian of F 

mag ; 3j function of 2d+l3 real variables on the same manifold. See equations ( l28l) 
and (129|) for explicit formulas. 

As before, the point Zj = fj, 7^ = 7^ is a critical point of the function /i™* ,^,^(2;). Similarly, 

Zj = fj, = is a critical point of the function /imag"°''(-2). 

For a fixed vector z in the vicinity of the eigenfunction / we find the critical points of functions 
^7ir--,7^j^^) and /imag""''(-2) with respect to the corresponding set of parameters. We choose the critical 
points close to (71, ... , 7^) and (0, . . . , 0) correspondingly. 

Concentrating first on /i™* (2;) we observe that now the terms that depend on a 7 have the form 

(37) T(7) = -7k„P-k.lV7- 

The critical point is g{z) = \zy\/\zu\ and it is a maximum or minimum depending whether 7 > or 
7 < (the sign of 7 is determined by the sign of the corresponding 7). Similarly to the proof of 
Theorem 15. 2[ the index of the critical points that are close to (71, . . . ,7/3) is /3 — p, where p is the 
number of negatives among 7^ . On the critical manifold the term (!37ll evaluates to 

(38) T((7) = -2sign(7)|z„||2;,|. 

Moving on to h'^ag'°'^{z) we first investigate its dependence on a single phase a. The terms 
involving a are of the form 

(39) Tmia) = -^e*"z2 - ^e"*"2;„ = -2|2;„||2;„| cos(a - 9uv), 
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where 9^^ = aTg{zu/zy). Since Zu,Zy G C are close to fu,fv ^ ^ correspondingly, the angle 9uv is 
close to if 7 = fy/ > and close to vr if 7 = f^/ /„ < 0. The corresponding critical point a = a{z) 
is then minimum (7 > 0) or maximum (7 < 0). 

Considering now the critical point of /imag""'^(2) as a function of all parameters ai, . . . ,a/3 (the 
variable z is fixed), we find that the Hessian is diagonal and therefore the index is the number of 
coordinate-wise maxima, that is p. Moreover, on the critical manifold the term ( l39l) takes the form 
identical to equation ( 138|) . Comparing the forms and /imag (2;) we conclude that, for 

each value of z, they coincide at their respective critical points, assuming the value 

h{z) = ^qu\zy,\'^ - ^ 2Re{z:;;Zy)- ^ 2sign{^j)\zu\\zy\. 

u {u,v)gE\S ej=(u,v)eS 

Thus indices of the point z = f with respect to the corresponding critical manifold of h'^^ -ypi^) 
and /imag '"''(^) are equal. This establishes a bridge between the two forms. Note that the bridging 
function h{z) is different from the form h{z) of the original graph F. The relations between the 
quadratic forms and the eigenvalues as functions of parameters are sketched in Fig. 15.31 We will 
now be repeatedly applying the reduction of critical manifolds. Theorem 14.51 The indices of the 
reductions are noted along the arrows in Fig. 15.31 (see below for explanations). 

Fixing the parameters 71, . . . , 7^, we find the critical point of h'^'^^ -yfji^) ^ function of z. Locally 
around the point zj = fj, 'jk = Ik, the critical point is the (0„ — p + l)-th eigenvector of the graph 
F™* According to Theorem 14.41 its index is 2(</)„ — p). We apply the reduction theorem to the 
result of Theorem 15.21 to conclude that the index of the point 2; = /, 7 = 7 with respect to all 
parameters is 

{n-l + l3-(f)n) + 2{(j)n -p) = n-l + (3 + (f)n-2p. 

Now we go down to the "bridge" critical point and then back up to the critical point z = f , a = 
of the function hmag'"'^ (z), applying the reduction theorem both ways. For the latter critical point 
we obtain the index 

{n-l + l3 + (j)n-2p)-{/3-p)+p = n-l + (f)n. 

Finally, for every choice of parameters ai, . . . , a^j, we find the critical point with respect to z that is 
close to z = f. This critical point /mag is the n-th eigenvector of the graph Fmag""'' (since / is the 
?7,-th eigenvector of F^J^j^g'"). According to Theorem 14.41 it has index 2(n— 1). We apply the reduction 

theorem one last time to conclude that the point (0, . . . , 0) as a critical point of /imag ''°''(/mag '"'') 
has index 

n - 1 + (pn -2{n - 1) = (pn - {n - 1). 
Since /imag ""''(/mag "'"") = A^i (Fmag ) , this coucludcs the proof of the theorem. □ 

6. Discussion 

Perhaps the most important feature of Theorem 1 1.1 1 is that it allows to access some of the features 
of the eigenfunction via the behavior of the corresponding eigenvalue under perturbation. It is known 
that the eigenvalues of the Laplacian are connected to the statistics of the closed paths on the graph. 
The connection is given through the so-called "trace formulae" , which can be obtained from a graph 
analogue of the Selberg zeta function, the Ihara zeta function [32l O HS]. An extension by Bartholdi 
[1] (see also [36]) was used in [39] to obtain a family of trace formulae including the ones for the 
magnetic Laplacian. Thus, the closed paths on the graph determine the spectrum of the magnetic 
Laplacian which, in turn, determines the nodal count. This, in principle, establishes the existence 
of a general connection between the nodal count and the closed paths. However, we are not aware 
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of any concrete general formulas. We note that such a connection has been earlier conjectured by 
Smilansky, with special cases reported in [241 [1] • 

We note that the eigenvalue A„(rmag"'"'') featuring in this paper is a well-studied object. It is the 
dispersion relation for the maximal Abelian cover of the graph F. One of the interesting questions 
regarding this object is the "full spectrum property" [29l |28l HI] : whether the continuous spectrum of 
the cover graph of a regular graph — in our terms, the union of ranges of the functions A„(rmag''"'') 
— contains no gaps. This question can be reformulated in terms of eigenfunctions of graphs rmag'""'' 
with all aj = or vr that have minimal and maximal number of sign changes. 

This, in turn, is related to the question of whether the extrema of the dispersion relation are 
always achieved at the symmetry points (namely, all aj = or vr). Examples to the contrary have 
been put forward in p6l [2T] . However, an important question remains, how can one characterize 
the extremal points that are not points of symmetry? In this direction, the duality with the cut 
graphs (Section [3]) might provide some answers. One can speculate that critical points of dispersion 
relation correspond to critical points of the eigenvalues of the cut graph that do not give rise 

to the eigenfunction of the graph F. Further, we conjecture that these "unclaimed" critical points 
correspond to eigenfunctions of F modified by enforcing Dirichlet conditions at some vertices. 

The results of the present paper are derived under the assumption that the eigenvalue is non- 
degenerate. While this is the generic situation with respect to the change in the potential Q, it is 
also interesting to consider what happens in the degenerate case. Linear Zeeman effect (the magnetic 
perturbation splits eigenvalues) suggests that the singularities of A„(Fmag are conical. It should 
be possible to define the index of the singularity point that does not rely on differentiability. 

Finally, it would be most interesting to generalize the results of the present paper to domains of 
W^. However, we immediately encounter a conceptual problem — the "number" of zeros is infinite. 
Still, some measure of instability of the eigenvalue under magnetic perturbation should be related to 
some measure of the zero set of the corresponding eigenfunction. This can be intuitively visualized 
by approximating the domain eigenfunction by eigenfunctions of a discrete mesh. 
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